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We report a new measurement of the ratio h/m-nb between the Planck constant and the mass 
of 87 Rb atom. A new value of the fine structure constant is deduced, oT 1 = 137.035 999 037(91) 
with a relative uncertainty of 6.6 x 10 -10 . Using this determination, we obtain a theoretical value 
of the electron anomaly a e = 0.001 159 652 181 13(84) which is in agreement with the experimental 
measurement of Gabrielse (a c = 0.001 159 652 180 73(28)). The comparison of these values provides 
the most stringent test of the QED. Moreover, the precision is large enough to verify for the first 
time the muonic and hadronic contributions to this anomaly. 
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The fine structure constant a characterizes the 
strength of the electromagnetic interaction. This dimen- 
sionless quantity is defined as: 



(1) 



where eo is the permittivity of vacuum, c the speed of 
light, e the electron charge and H the reduced Planck 
constant (h = h/2ir). It appears in the expressions of 
the ionization energy of hydrogen atom, of the fine and 
hyperfine structures of atomic energy levels, and it is the 
parameter of the quantum electrodynamics (QED) calcu- 
lations. Its measurement in different domains of physics 
is a test of the consistency of the theory. The most accu- 
rate value is deduced from the combination of the mea- 
surement of the electron anomaly a c with a very difficult 
QED calculation. The last result, by Gabrielse at Har- 
vard University, gives a value of a with a relative uncer- 
tainty of 3.7 x 10" 10 [H, Q. Nevertheless this impressive 
result is fully dependent on QED calculations and can be 
liable to a possible error: in 2007, Aoyama et al detected 
an error which shifted the a value by 4.7 x 10~ 9 0-01 • 
Consequently, to check these calculations, another deter- 
mination of a is required. Up to now the other QED 
independent determinations of a were less accurate by 
at least an order of magnitude. The measurement of the 
quantum Hall effect provides an a value with an uncer- 
tainty of 1.8 x 10~ 8 |5J] and the accuracies of the deter- 
minations deduced from the recoil velocity measurement 
were respectively 7.7 x 10~ 9 and 4.6 x 10~ 9 for the Cesium 
and Rubidium experiments 0, 0] • 

In this Letter we present a new measurement of the 
ratio /i/mRb between the Planck constant and the mass 
of 87 Rb atom and we obtain a new value of a: 



a" 1 = 137.035 999 037(91). 
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With a relative uncertainty of 6.6 x 10 10 , this value im- 
proves our precedent result by a factor of about seven Q ■ 
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FIG. 1: Determinations of a with a relative uncertainty 
smaller than 10 -8 ; a c (UW): measurement by Dehmelt at the 
University of Washington @|; h/mc B - measurement of the 
Cesium recoil velocity at Stanford h/mnb' measurement 
of the Rubidium recoil velocity at Paris in 2006 [f| and 2008 
0; o e (Harvard): measurement of g — 2 at Harvard University 
in 2006 [3 and 2008 [J; CODATA 2006: best adjustment by 
the Committee on Data for Science and Technology 0; the 
arrow corresponds to the shift of the values of a due to the 
reevaluation of the QED calculation of a c in 2007 @, H|. 



The comparison with the value deduced from the elec- 
tron anomaly provides the most stringent test of QED 
[l|. Indeed there is a very good agreement with this 
last value (o^ 1 = 137.035 999 084(51)) as illustrated on 
Fig. [TJ This agreement confirms together the recent g — 2 
measurement of Gabrielse by comparison with the value 
obtained by Dehmelt at the University of Washington Q 
and the recent correction found in the calculation of the 
electron anomaly Q. The discussion on this agreement 
will be presented at the end of this Letter. 

The fine structure constant is deduced from the mea- 
surement of /i/mRb thanks to the relation: 
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where m c is the electron mass. In equation ([3]), the 
Rydbcrg constant Roc and the mass ratio mnb/wje are 



15567824.42 +/- 0.15 Hz 



-15567822.07 +/- 0.16 Hz 



known with an accuracy of 7 x lCr 12 H M, till and 
4.4 x 10~ 10 [ill, 13| respectively: the limiting factor is 
the ratio h/m^b- In our experiment, it is deduced from 
the measurement of the recoil velocity v r = hk/m^b of a 
Rb atom when it absorbs a photon of momentum hk. 

The principle of the experiment has been described 
previously Ml. It combines a Ramscy-Borde atom inter- 
ferometer [14| with Bloch oscillations (BO). The method 
is to coherently transfer many recoils to the atoms at rest 
and to measure the final velocity distribution. The ex- 
periment develops in three steps, i) Firstly, a pair of 7r/2 
pulses of a Raman transition transfers the 87 Rb atoms 
from the F = 2 hypcrfine sublevel to the F = 1 one and 
produces a fringe pattern in the velocity distribution of 
these atoms. The width of the envelope of this velocity 
distribution varies inversely with the 7r/2 pulse duration 
r, while the fringe width varies as 1/Tr, where Tr, is 
the delay between the two 7r/2 pulses, ii) Secondly, we 
transfer to the selected atoms as many recoils as possi- 
ble by means of BO. Bloch oscillations have been first 
observed in atomic physics by the groups of Salomon in 
Paris and Raizen in Austin 15l-ll7|. They can be inter- 



preted as Raman transitions in which the atom begins 
and ends in the same energy level, so that its internal 
state (F = 1) is unchanged while its velocity has in- 
creased by 2v r per BO. Bloch oscillations are produced 
in a one dimension vertical optical lattice which is ac- 
celerated by linearly sweeping the relative frequencies of 
the two counter propagating laser beams (frequencies V\ 
and v 2 )- This leads to a succession of rapid adiabatic 
passages between momentum states differing by 2hk. iii) 
Finally, the final velocity of the atoms is measured by a 
second pair of 7r/2 pulses which transfers the atoms from 
the F = 1 to the F = 2 hyperfinc level. The frequency 
difference between the two pairs of 7r/2 pulses is scanned 
to obtain a fringe pattern from which the velocity varia- 
tion between the two pairs of 7r/2 pulses is deduced. 

The experimental setup uses a double vacuum cell. A 
two dimensional magneto optical trap (2D-MOT) pro- 
duces a slow atomic beam (about 10 9 atoms/s at a ve- 
locity of 20 m/s) which loads during 250 ms a three di- 
mensional magneto optical trap. Then a a + — o~ mo- 
lasses generates a cloud of about 2 x 10 s atoms (in the 
F = 2 hypcrfine level) with a 1.7 mm radius and at a 
temperature of 4 /xK. A vertical magnetic field of 7 /jT 
is applied and a radio frequency pulse is used to select 
atoms in the F = 2, mp = state: a first pulse transfers 
the atoms from the F = 2, mp = Zeeman sub level to 
the F = 1, mp = one. A laser beam pushes away the 
atoms left in the F = 2 hyperfine level and the atoms 
in the F = 1, mp = level come back to the F = 2, 
m p = level with a second radio frequency pulse. Then 
we follow the procedure described in reference 0. The 
atoms are accelerated and decelerated by two sequences 
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FIG. 2: Example of the four spectra needed to deduce 
h/rriYib (see the text). Each spectrum represents the quan- 
tity N2/(Ni + N2) where Ni and N2 are the population of 
the F = 1 and 2 levels in function of the frequency difference 
between the two pairs of n/2 pulses. The measured position 
of the central fringe is indicated above the spectra. 



of 300 BO (duration of 4.6 ms with a delay between them 
of 10.3 ms) to make an atomic elevator and displace the 
atomic cloud in the upward or downward direction. The 
atoms are accelerated in the opposite direction with 500 
BO (duration 5.6 ms) before the first pair of 7r/2 pulses 
which produces the Raman transition from the F = 2, 
mp = level to the F = 1, mp = one (delay Tr = 
10 ms and duration r = 600 jits). The atoms left in the 
F = 2 are pushed away. Finally, 500 BO decelerate the 
atoms before the second pair of 7r/2 pulses and the pop- 
ulation of the F = 2 and F = 1 levels are measured 
with a time of flight technique. The laser beams used 
for the Raman transitions and the Bloch oscillations are 
sent on the atoms thanks to two optical fibers following 
the scheme described in reference [9| . They are very well 
collimated with a beam waist of 3.6 mm and the power 
of each beam used for the BO is about 150 mW. 

A value of h/m^b is obtained by recording four spec- 
tra. To cancel the velocity variation gT due to the grav- 
ity g (T = 19 ms is the delay between the two pairs 
of 7r/2 pulses), the atoms are accelerated alternatively 
upward and downward and the difference between the 
results eliminate gT. Moreover, for each initial accel- 
eration, two spectra are recorded by exchanging the di- 
rections of the Raman beams to eliminate the parasitic 
level shifts due to the Zeeman effect or to the light shifts. 
The figure shows an example of records. The acqui- 
sition time is about five minutes. For each spectrum, 
the Doppler shift is obtained with a relative accuracy of 
about 10~ 8 . The ratio h/m^b can be then determined 
from: 
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FIG. 3: Measurements of the ratio h/m-Rb during about 15 
hours. The standard deviation of the mean is 4.4 x 10~ 10 
with \ 2 /( n ~ 1) = 1-05. The inset shows the autocorrelation 
function of these 170 measurements. The solid and dashed 
lines represent the 1 a and 2 a standard deviation of the 
autocorrelation function. 



TABLE I: Error budget on the determination of 1/q (system- 
atic effect and relative uncertainty in part per 10 10 . 



Source 



Correction 



Relative 
uncertain ty 



Laser frequencies 




1.3 


Beams alignment 


-3.3 


3.3 


Wavefront curvature and Gouy phase 


-25.1 


3.0 


2nd order Zeeman effect 


4.0 


3.0 


Gravity gradient 


-2.0 


0.2 


Light shift (one photon transition) 




0.1 


Light shift (two photon transition) 




0.01 


Light shift (Bloch oscillation) 




0.5 


Index of refraction atomic cloud 






and atom interactions 




2.0 


Global systematic effects 


-26.4 


5.9 


Statistical uncertainty 




2.0 


Rydberg constant and mass ratio 




2.2 


Total uncertainty 




6.6 



where N = 500 is the number of Bloch oscillations in 
both opposite directions, ks is the Bloch wavevector and 
ki and ki are the wavevectors of the two Raman beams. 
Consequently, the spectra of Fig. [5] give h/m^b with a 
relative statistical uncertainty of 5 x 10~ 9 (2.5 x 10~ 9 
for a). We have recorded about 1370 spectra analog 
to the figure [H An example of 170 determinations of 
h/m^b obtained during one night is displayed on Fig. [3] 
The autocorrelation function of these measurements (see 
the reference (HI), which is reported in the inset of Fig. 
[H shows no correlation between the successive measure- 
ments of h/mnt,. 

Table U gives the error budget. The systematic ef- 
fects are reduced compared to our previous measure- 
ments 0, Q. The lasers arc locked on a Fabry- Perot 
cavity stabilized with a standard laser and their frequen- 
cies are frequently measured with a frequency comb to 
reduce the frequency uncertainties at less than 50 kHz. 
The maximum angle between the lasers used for the Ra- 
man transitions and the Bloch oscillations is estimated 
to 40 fir&d from the coupling between the two optical 
fibers. Moreover this value has been confirmed by the 
observation of the effect of the misalignment between the 
Bloch beams (see Fig. H]). The effect of the Gouy phase 
and the wave front curvature, which has been reduced 
by increasing the waist of the laser beam from w = 2 
mm to w — 3.6 mm, has been carefully controlled with 
a Shack-Hartmann wave front analyzer. The parasitic 
magnetic field has been reduced with a double magnetic 
shield and a precise mapping of the magnetic field gives 
now a relative correction of 4 x 10~ 10 . Thanks to the 
good collimation of the laser beams, the section of the 
laser beams varies by about 4 x 10~ 3 along the atomic 
trajectory and the result is a very good cancelation of the 
light shift effects between the upward and downward tra- 
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FIG. 4: (Color online). Sensitivity of the /i/rngb measure- 
ments to the alignment: results of the measurements when 
the angle between the two Bloch beams is modified. The 
most precise point corresponds to the measurements in Fig. 
[3] The difference between this value and the summit of the 
parabola which is fitted to the data is 0.35 ppb. 

jectory. From the density of the cloud of cold atoms after 
the RF selection and the two first sequences of BO (about 
2 x 10 s atoms/cm 3 ), the effects of the refractive index and 
of the interactions between the atoms are estimated at a 
10~ 10 level, corresponding to a conservative uncertainty 
of 2 x lO" 10 in Table U Thanks to the double cell with 
a differential pumping, the effect of the refractive index 
due to the background vapor (about 10 7 atoms/cm 3 ) is 
now at the negligible level of few 10 -11 . 

Taking into account all these corrections, the measured 
value of h/m Rh is 4.591 359 2729 (57) x lO^nrV 1 , and, 
using the values of Roc , m c [H and mRb [IH , we obtain 
the value of a given by equation p|). 

With this new result the QED can be tested at a level 
better than 10~ 9 . During two decades, the theory of 
a c has been improved by Kinoshita and collaborators 0, 
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provides the most stringent test of QED. In the future, 
the sensibility of the interferometer can be increased by 
using a larger area interferometer 22 , 23[ and the experi- 
ment can be improved to reduce the correction due to the 
Gouy phase and the wave front curvature and divide the 
uncertainty by a factor of two. Then a main limitation 
will be the uncertainty of the mass ratio m^b/we- 

This experiment is supported in part by IFRAF (Insti- 
tut Francilien de Recherches sur les Atomes Froids), and 
by the Agence Nationale pour la Recherche, FISCOM 
Project- ( ANR-06-BL AN-0 1 92) . 



FIG. 5: (Color online). Upper figure: in blue, relative con- 
tributions to the electron anomaly of the different terms of 
equation ((5]), in red their uncertainties. The dashed line 
corresponds to the relative uncertainty of the new a value. 
Lower figure: comparison of the measurements of the elec- 
tron anomaly (UW 1987 @] and Harvard 2008 Q) with the 
theoretical value obtained by using the new value of a (la- 
bel Rb 2010). The point "Rb 2010 - only QED" is obtained 
without the last term of equation ([5j|. 



Il9l . |2C§. The anomaly is expressed as a sum of terms 
in power of a/n and of additive terms which take into 
account the contributions due to the muon, the tau, the 
weak interaction and the hadrons: 
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This equation and the new value of a give a theoreti- 
cal value of the anomaly a c = 0.001 159 652 181 13(84) 
where the uncertainty is due to the theory (33 x 10 -14 ) 
and to the a measurement (78 x 10 -14 ). The compari- 
son with the experimental result [H gives the difference 
ac(thcory) - a c (cxp.) = (40 ± 89) x 10~ 14 . The rela- 
tive agreement between the experiment and the theory 
is at the level of 7.7 x 10~ 10 . The figure [5] shows the 
contributions of the different terms of equation (up- 
per part) and the comparison between the experimen- 
tal and theoretical values (lower part). The accuracy of 
the a measurement is sufficient to test for the first time 
the contributions due to the muon and hadrons. If we 
suppose the exactness of the QED calculation, this very 
good agreement provides a strong limitation to a possible 
structure of the electron [2(| or to the existence of new 
dark matter particles (ilj . 

In conclusion we have presented a recoil- velocity mea- 
surement of Rubidium and obtained a new determination 
of the fine structure constant with a relative uncertainty 
of 6.6 x 10~ 10 . The combination of this result with the 
measurement and the calculation of the electron anomaly 



[1] D. Hanneke, S. Fogwell and G. Gabrielse, Phys. Rev. Lett. 

100, 120801 (2008). 
[2] T. Aoyama, M. Hayakawa, T. Kinoshita and M. Nio, 

Phys. Rev. Lett. 99, 110406 (2007). 
[3] B. Odom, D. Hanneke, B. D'Urso and G. Gabrielse, Phys. 

Rev. Lett. 97, 030801 (2006). 
[4] G. Gabrielse, D. Hanneke, T. Kinoshita, M. Nio and 

B. Odom, Phys. Rev. Lett. 99, 039902 (2007). 

[5] P. J. Mohr, B.N. Taylor and D.B. Newell, Rev. Mod. 

Phys. 80, 633 (2008). 
[6] A. Wicht, J.M. Hensley, E. Sarajilic and S. Chu, Physica 

Scripta T102, 82 (2002). 
[7] M. Cadoret, E. de Mirandes, P. Clade, S. Guellati- 

Khelifa, C. Schwob, F. Nez, L. Julien, and F. Biraben, 

Phys. Rev. Lett. 101,230801 (2008). 
[8] R.S. Van Dyck, P.B. Schwinberg and H.G. Dehmelt, 

Phys. Rev. Lett. 59, 26 (1987). 
[9] P. Clade, E. de Mirandes, M. Cadoret, S. Guellati- 

Khelifa, C. Schwob, F. Nez, L. Julien and F. Biraben, 

Phys. Rev. A 74, 052109 (2006). 
[10] Th. Udem, A. Huber, B. Gross, J. Reichert, M. 

Prevedelli, M. Weitz, and T. W. Hansch, Phys. Rev. Lett. 

79, 2646 (1997). 
[11] B. de Beauvoir, C. Schwob, O. Acef, J. -J. Zondy, L. Joze- 

fowski, L. Hilico, F. Nez, L. Julien, A. Clairon and 

F. Biraben, Eur. Phys. J. D. 12, 61 (2000). 
[12] MP. Bradley, J.V. Porto, S. Rainville, J.K. Thompson, 

and D.E. Pritchard, Phys. Rev. Lett. 83, 4510 (1999). 
[13] B.J. Mount, M. Redshaw, and Edmund G. Myers, Phys. 

Rev. A 82, 042513 (2010). 
[14] Ch.J. Borde, Phys. Lett. A 140, 10 (1989). 
[15] M. Ben Dahan, E. Peik, J. Reichel, Y. Castin and C. Sa- 
lomon, Phys. Rev. Lett. 76, 4508 (1996). 
[16] E. Peik, M. Ben Dahan, I. Bouchoule, Y. Castin and 

C. Salomon, Phys. Rev. A 55, 2989 (1997). 

[17] S. R. Wilkinson, C. F. Bharucha, K. W. Madison, Qian 
Niu, and M. G. Raizen, Phys. Rev. Lett. 76, 4512 (1996). 

[18] T.J. Witt, Metrologia 44, 201 (2007). 

[19] T. Kinoshita and M. Nio, Phys. Rev. D 73, 013003 
(2006). 

[20] G. Gabrielse, D. Hanneke, T. Kinoshita, M. Nio, and B. 

Odom, Phys. Rev. Lett. 97, 030802 (2006). 
[21] C. Boehm and J. Silk, Phys. Lett. B 661, 287 (2008). 
[22] S.-W. Chiow, S. Herrmann, S. Chu, and H. Miiller, Phys. 

Rev. Lett. 103, 050402 (2009). 
[23] P. Clade, S. Guellati-Khelifa, F. Nez, and F. Biraben, 

Phys. Rev. Lett. 102 240402 (2009). 



